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ABSTRACT 

AGN heating, through massive subrelativistic outflows, might be the key to solve 
the long-lasting 'cooling flow problem' in cosmological systems. In a previous paper, 
we showed that cold accretion feedback and, to a lesser degree, Bondi self-regulated 
models are in fact able to quench cooling rates for several Gyr, at the same time 
preserving the main cool core features, like observed density and temperature profiles. 
Is it true also for lighter systems, such as galaxy groups? The answer is globally yes, 
although with remarkable differences. Adopting a modified version of the AMR code 
FLASH 3.2, we found that successful 3D simulations with cold and Bondi models are 
almost convergent in the galaxy group environment, with mechanical efficiencies in the 
range 5 x 10"'* — 10"'^ and 5 x 10"^ — 10~^, respectively. The evolutionary storyline of 
galaxy groups is dominated by a quasi-continuous gentle injection with sub-Eddington 
outflows (with mechanical power P ^ 10''* erg s~*, v ^ 10* km s^*). The cold and 
hybrid accretion models present, in addition, very short quiescence periods, followed by 
moderate outbursts (10 times the previous phase), which generate a series of 10—20 kpc 
size cavities with high density contrast, temperatures similar to the ambient medium 
and cold rims. After shock heating, a phase of turbulence promotes gas mixing and 
diffusion of metals, which peak along jet-axis (up to 40 kpc) during active phases. At 
this stage the tunnel, produced by the enduring outflow (hard to detect in the mock 
X-ray surface brightness maps), is easily fragmented, producing tiny buoyant bubbles, 
typically a few kpc in size. In contrast to galaxy clusters, the AGN self-regulated 
feedback has to be persistent, with a 'delicate touch', rather than rare and explosive 
strokes. This evolutionary difference dictates that galaxy groups are not scaled-down 
versions of clusters: AGN heating might operate in different regimes, contributing to 
the self-similarity breaking observed. 

Key words: cooling flows - galaxies: active - galaxies: jets - hydrodynamics - 
intergalactic medium ~ X-rays: galaxies: groups. 



1 INTRODUCTION 

Remarkably, most of the galaxies in the universe reside in 
small groups (up to 70% in the nearby zone; Geller & Hu- 
cra||1983 Mulchaey '20001 . Nevertheless, relatively little at- 



tention has been devoted to groups with respect to their 
'big brothers', galaxy clusters. Groups, with typical masses 
around 10^'^Mq, are thus a key class of objects, containing 
a significant fraction of the overall universal baryon budget 
and being a nursery for massive (early-type) galaxies. 

Galaxy groups are X-ray sources, with diffuse extended 
(50-500 kpc) emission similar to that of rich clusters, but 
with lower luminosities, ranging from 10'*^ erg s~^ up to 
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several times 10*^ erg s ^ , with temperatures around ~ 0.3- 
2 keV ( |Mulchaey|[2550l ). 



It is known that clusters and groups do not follow 
the global scaling relations predicted by self-similar mod- 
els (e.g. |Kaiser|p86l [Evrard fc Henr^poT] [Osmond fc| 



Ponman 



2004 Vikhlinin et al. 20061. Groups are not per 



feet scaled down versions of clusters. If groups are sampled 
at sufficiently large radii, some X-ray scaling relations (e.g. 
luminosity— temperature, mass— temperature) seem similar 



to those followed by massive clusters ( jOsmond fc Ponman 
2004[|Sun et al.|2009[ ), with the scatter increasing for T < 1 
keV. However, systematic differences between groups and 
clusters exist, with the former having clearly lower gas (and 
baryon) fractions and flatter entropy profiles ( [Finoguenov et| 
al.|2005|[Voit|2005((Sun et al.|2009[ ). Moreover, groups have 
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velocity dispersions of several hundreds km s ^ , compara- al.|2001 2003 Tamura et al. 2003 Peterson fc Fabian|2006 1 



ble to the internal velocity dispersion of individual galaxies. 
For this reason, some processes are more (merging) or less 
(ram-pressure stripping, galaxy harassment) relevant than 
in clusters. 

The discrepancies between the observed scaling rela- 
tions and the predictions of self-similar models, and espe- 
cially the systematic differences between groups and clus- 
ters quoted above, can be explained by the presence of 
non-gravitational heating, acting during or after the gravi- 
tational collapse. In lower mass systems the same amount of 
energy per unit mass has a stronger effect. This is a key point 
for the mechanical AGN feedback problem investigated here 



and in our previous work (Gaspari et al. 2011 Gil there- 
after) . 

As in clusters and elliptical galaxies, many galaxy 
groups display the typical signatures of cooling flowQ (a 
most notable example being NGC 5044). In the center, 
where the radiative cooling time is the shortest (~ 10^ yr), 
the temperature of the intragroup medium (IGM) drops 
by nearly 50%. Moreover, the X-ray surface brightness is 
sharply peaked, with the maximum coincident with a mas- 
sive elliptical galaxy. 

A large sample of X-ray groups is presented in |Sun et| 
al. (20091. Gathering 43 objects with Chandra data, they 



confirmed the existence of a 'universal' temperature profile 
( De Grandi fc Molendi[2002{ [VoitpOOS} [Leccardi fc Molendi 



20081). The group profiles are slightly steeper than those of 
clusters for r ^ 0.15 r^oo- In the inner cool core, the scatter 
is significant: the majority presents positive gradients, while 
few are approaching a fiat temperature profile (a possible 
signature of a heating cycle). Sun et al. produced also the 
entropy profiles, which retain the thermal history of the gas. 



They again confirmed previous important works (Ponman 
et al.|199 9, 2003 : 'Fino guenov et al.|2002[ ), suggesting devia- 
tions in the entropy from the self-similar relation {K oc T), 
showing a clear fiat core, and an excess relative to massive 
clusters. 

One of the most striking riddles of galaxy group astro- 
physics concerns the so called 'cooling flow problem'. This 
paper attempts to solve it, finding plausible consistent mod- 
els, as has been done for galaxy clusters in Gil. The classical 
cooling fiow (CF) model ( Fabian|1994 for a review) predicts 
that, as the gas radiates, its entropy decreases, and thus gets 
compressed by the surrounding gas, causing it to subsoni- 
cally flow inward. The estimated (isobaric) cooling rates is 
given by: 

2 fimp 



M 



5 kBT 



(1) 



where Lx is the X-ray luminosity, /j, the mean molecular 
weight, fee the Boltzmann constant, rrip the proton mass. 
For typical luminosities of groups (~ 10^'' erg s~^), the 
relation implies cooling rates of several tens M0/yr. Such 
high predicted M have been in contrast for years with the 
estimates of star formation rates and cold gas masses (at 
least an order of magnitude lower) , but the definitive obser- 
vational evidence against the traditional CF model has been 



the spectral analysis of XMM-RGS spectra (e.g. Peterson et 



Interestingly, the analogues of non-cool core clusters on the 
group scale have been rarely detected i Johnson et al.|2011 1 . 



In the last decade the new generation XMM- 
Newton and Chandra telescopes have radically changed the 
panorama of the evolution of hot gas in groups and clusters. 
ROSAT-HRI and, especially, Chandra high resolution X-ray 
images show clear evidence of AGN-gas interaction from gi- 



Blanton et al.||2001| |Finoguenov & 


: Jones||2001| |Jones et al.| 


2002 


McNamara fc Nulsen 


2007 


and references therein). 



Evidences of such activity have also been identified in many 
groups ([Allen et al . "20(36: M orita et al]|2006| |Jetha et ah] 



2008 Gastaldello et al.,,2009j . The fairly common presence 
of X-ray cavities, often coincident with lobes of radio emis- 
sion connected to the core of the central galaxy by a radio 
jet, indicates that AGN inject energy in the IGM in kinetic 
form (outflows) and as relativistic particles. 

Moreover, from dozens of giant ellipticals (gE) () Nulsen 



et al.|2007| ) to rich clusters ( [Rafferty et aLp 006 2008| there 



is a continuous scaling relation between the AGN power, 
associated with the cavity formation, and the core X-ray lu- 
minosity. This trend is more evident in low mass systems. 
Although this connection is likely to be a necessary but not 
a sufficient requirement for a successful heating scenario, it 
strongly suggests that the dominant heating process man- 
ifests itself generating bubbles in the ICM, through some 
kind of 'directional' energy input, such as jets or coUimated 
outflows (rather than spherically symmetric forms). 

AGN outflows seem to be a very promising mechanism 
to simultaneously explain the quenching of cooling flows 
(and star formation), the deviation of scaling relations from 
self-similar laws and the presence of X-ray cavities and 
weak shocks in the intergalactic (interstellar) medium. As 
anticipated before, the impact of this type of heating is 
much more evident in less bounded systems, such as groups. 
We can therefore speculate that gentler outflows are more 
appropriate to regulate the thermodynamic evolution of the 
IGM, while strong bursts should easily deform the main 
properties. 

Regarding clusters, we can account for several varie- 



gated computations (e.g. Ruszkowski et al.|2004[ Omma et 
'al. 2004 ; Bruggen et al."2005' 'Brighonti & Matliews"2006| 
^Sternberg et al.,2007( |Gaspari et al.^2009, 2011 among oth- 
ers), but little work has been done to explore the role of 
AGN heating on the intragroup medium, especially from a 
theoretical and numerical point of view. [Brighe nti fc Math-| 
( 20 035 studied thermal heating (plus conduction) in 2D 



simulations, showing that this type of feedback, while very 
efficient in stopping the cooling process, often generates neg- 
ative temperature gradients, with the total erasure of galac- 
tic cooling fiows commonly observed in gE ( [Sun et al.|2005 



200 7|. C arrying out cosmological simulations ( [Puchwein et| 
|al.|2008 McCarthy et al. 2010 1 does not permit, nowadays, a 
detailed analysis of the feedback dynamics in the inner cool 
core, mainly because of the lack of resolution, and are thus 
more suited to study global sample-averaged properties. Fur- 
ther, none of the above mentioned works has investigated a 
purely mechanical feedback. 

In this paper we will therefore describe, for the first 
time, 3D self-regulated hydrodynamical simulations of AGN 
outflows in an exemplary galaxy group, lasting at least 7 
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Gyr. The positive results of similar simulations for rich 
galaxy clusters (Gil) are a strong motivation for a detailed 
analysis of this (mechanical) feedback mechanism in lower 
meiss systems. 

As in Gil, we assume that AGN outflows are the key 
element in solving the cooling flow problem. The physical 
foundation is that the relativistic jet, produced by the active 
nucleus, entrains ambient gas, strongly decelerating within 
a scale-length of just a few kpc. Such subrelativistic AGN 
outflows have been observed in several cD galaxies in every 
band, thanks to blueshifted absorption lines: optical (|Nes- 
vadba et aL]|2008[ [20lT] ), U V and X-Ray (see |Crenshaw et" 



al.|2003| fora review), 21-cm ( [Morganti et al.|2605[|2007 



also the references in Gil). They occur at a distance of a 
kpc with velocities around 10^ — 10* km s"'^. The geometry 
of the bipolar outflows is still unclear, with usually narrow 
opening angles, and possibly some kind of precession. 

The timing of the feedback is also an important issue. 
From observations and from a 'thermostat' assumption the 
duty cycle has to be similar to the cooling time, around 10^ 
yr. We will simplify the feedback self-regulation, testing hot 
(Bondi) or cold accretion models with different mechanical 
efflciencies. Fully 3D simulations are required to compute 
the chaotic and turbulent flow, with associated instabilities. 

In summary, this work will focus on a set of heated cool- 
ing flow models with a variety of AGN feedback mechanisms, 
aimed to answer the following key question: are AGN out- 
flows able to prevent the IGM from cooling and at the same 
time preserve the observed cool core appearance for many 
Gyr? 



2 THE COMPUTATIONAL PROCEDURE 

In order to study the behaviour of AGN outflows in galaxy 
groups, we carried out several simulations with a substan- 
tially modifled version of FLASH 3.2 ( |FryxeU et al.|2000| ), a 
popular 3D adaptive mesh refinement (AMR) code. FLASH 
uses the Message-Passing Interface (MPI) library to achieve 
portability and efficient scalability on a variety of different 
parallel high-performance computing systems. The simula- 
tions were usually run on 128 processors of IBM P575 Power 
6 (SP6) at CINECA supercomputing centre. 

We tested different numerical schemes to solve the hy- 
drodynamic Euler equations. As in Gil, we widely used 



the Piecewise-Parabolic Method (PPM) solver (Colella & 
Woodward|1984 1 , particularly appropriate to describe shock 
fronts. It uses directional splitting to advect quantities 



(Strang 19681. This makes the instantaneous injection of 



the outfiow not trivial, since we must manage two equal 
timesteps with a sweep x — y — z and then z — y — x. Thus, we 
decided to also adopt other unsplit methods (single timestep 
evolution). We tested Roe and Lax-Friedrichs, with different 
slope limiters (van Leer, minmod, etc.). They work pretty 
well for smooth fiows, but not for a powerful jet ignition, 
producing unstable results also with low Courant-Friedrichs- 
Lewy (CFL) numbers (< 0.1). On the contrary, the HLLC 
solver ( |Toro||1999 1 is very stable also for higher CFL num- 
bers (~ 0.5), with a good description of shocks and contact 
discontinuities. In the end, results with both methods are 
comparable, with PPM being more accurate (3rd order), 



while HLLC faster and more manageable. In many cases we 
preferred accuracy over speed. 

The innovations implemented in FLASH reside, mainly, 
in several source and sink terms added to the usual hydro- 
equations. In conservative form: 

dp 



V ■ [pv) = ap, - q- h S'ljet , 



dpv 



+ V ■ {pv(g,v) + VP = pgoM + S2,iot , 



dps „ r/ ^ I '"^ 

^ + V • [{pe + P)v]^ pv- + ap^ I ^" + T 
— noniA(r, Z) + Ssjct , 



P = (7-l)p(e-y 



(2) 
(3) 

(4) 
(5) 



where p is the gas density, v the velocity, e the specific total 
energy (internal and kinetic), P the pressure, the grav- 
itational acceleration (due to the dark matter halo plus the 
central elliptical, see below), and 7 = 5/3 the adiabatic in- 
dex. The temperature is computed from P and p using 
with an atomic weight p, ~ 0.62, appropriate for a totally 
ionized plasma with 25% He in mass. 

In Gil we describe in-depth all the source and sink 
terms (right hand side of the first three equations). Here we 
give a summary. 

Radiative cooling is modeled according to |Sutherland| 



& Dopita ( 1993| ) cooling function A(r, Z) for a fully ionized 



plasma (rie and ni are the number density of electrons and 
ions). For typical groups, the central metallicity Z is about 
solar (e.g., Rasmussen & Ponman 20091, with a negative 



gradient at larger radii. For simplicity, we set the abundance 
equal to 1 Zq. 

Supernovae la (SNIa) and stellar winds (SW), in the 
central elliptical galaxy are implemented following | BrighentTl 
[fc Mathe ws (2002), with a stellar mass loss rate a{t) domi- 
nated by stellar winds oc t~^'^. 

The key prescription for treating the cold gas is the 
dropout term (Eq. ([2|), proportional to q = 2exp(— (T/5 x 
10^)^). It removes the gas below T = 10* ^ K, whose physi- 
cal evolution cannot be properly followed by our cod«Q Be- 
cause it has been proved that the total mass of cooled gas is 
insensitive to the presence of the dropout term or its func- 



tional form (see Brighenti & Mathews 20001, we can still 
calculate the gas cooling rate accurately. 

Outfiow source terms 5i,2,3,jet will be explained for ev- 
ery type of feedback in Section 2.2. Injection will be done 
with two different methods: directly into the domain (with- 
out a mass infiow) or through boundary condition at 2: = 
(Si >0). 

We calculate the fiow evolution for 7 Gyr. Local cos- 
mological systems have very different formation times and 



are subject to relatively frequent mergers (Cohn & White 
|2005j ). We conservatively decided to run our simulations for 
such a long time to take into account the systems which do 
not undergo major mergers since z ~ 0.7. Moreover, cool 
cores in massive objects may survive major mergers (e.g., 

^ Tests without the dropout show the artificial assembly of cold 
clumps near the center. 
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Table 1. Parameters and properties of the most relevant models. 



while observations suggest the coohng rate is at least an 



order of magnitude lower, less than a few M0/yr (David 



Model 


Feedback 


efficiency (e) 


Notes 


CF 
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continuous 
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continuous 
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Bondi 
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cold timing 
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cold timing 


Bi5em2 
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5 


X 10-2 


cold timing 
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Bondi 




10-1 


cold timing 


C5eni5 


cold 


5 


X 10 ^ 
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cold 




10 
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cold 


5 


X 10-* 


. 


ClemS 


cold 




10-3 




IntSlOl 


5-10 Myr cycle 






Pj = 10-5PEdd 


IntSlOm 


5-10 Myr cycle 






Pj = lO-^PEdd 


IntSlOh 


5-10 Myr cycle 






Pj = 10-3PEdd 


Eth50 


thermal (-1- Be) 


5 


X 10-3 


50% Eth 


I01-9 


InOut (-1- Be) 




10-3 


Mln = O.lMacc 


IOen40 


InOut (-1- Be) 


5 


X 10-3 


40Mout (I03-7) 


lOenSO 


InOut (-1- Be) 


5 


X 10-3 


SOMout (I03-7) 



Burns et al.|2008 1 , implying that the central AGN feedback 
is continuously needed. 

Note that several feedback schemes could delay exces- 
sive gas cooling for only 1 Gyr before failing (see Gil). Thus, 
we stress that a few hundreds Myr evolution is not sufficient 
to make predictions on the heating versus cooling regulation 
and might lead to misleading extrapolations. 

2.1 The group model and initial conditions 

We choose NGC 5044 as the template for a typical X-ray 
bright group of galaxies: all the results we present should 
be relevant for any object in this category. NGC 5044 was 
one of the first cool core groups observed, due to its high 
X-ray brightness, with low redshift 0.009 (see Buote et al. 
[2003i [Gastaldello et al.1[2009l [David et al.||2009[ for recent 
Chandra and XMM data). The estimated X-ray bolomet- 
ric (0.1 — 100 keV) luminosity is around 2 x 10*3 erg s-^, 
with a virial mass of ~ 4 x 10'^3 Mq. The large scale mor- 
phology is very smooth and nearly spherical, despite a little 
disturbance in the form of a south-eastern cold front. At the 
contrary, the core (~ 10 kpc) has been strongly perturbed by 
recent outbursts from the central AGN. The group presents 
many small radio quiet cavities with a nearly isotropic dis- 
tribution and moderate mechanical power (~ 10*^ erg s-^). 
The two biggest cavities at larger radii (10 — 20 kpc) are 



instead filled with radio emission at 235 MHz (David et al 



2009 1 and the associated mechanical power seems to be able 
to balance radiative losses. The NGC 5044 X-ray image also 
shows several cold filaments coincident with Ha and dust 
emission, indicating a physical connection between the var- 
ious gas phases ( [Gastaldello et al.||2009| ). The GMRT ob- 
servation (610 MHz) reveals the presence of extended ra- 
dio emission with a torus-like morphology, threaded by the 
largest filament (probably cold material being uplifted from 
the center). 

The cooling time of the hot gas within the central 2 
kpc is just ~ 4 X 10^ yr. In the absence of feedback, a pure 
cooling flow model predicts tens of M©/yr (see Eq. (1)), 



et al.|[2009| ). This is the key requirement for our computed 



models 

This galaxy group consists of a luminous giant elliptical 
galaxy (NGC 5044) surrounded by a cluster of ~ 150 low 
luminosity dwarf galaxies, mostly of early type. Thus, we 
choose to model the elliptical galaxy with a de Vaucouleurs 
profile ( [Mellier fc Mathez|[1987[ ) with total stellar mass M<, 
~ 3.4 X 10 Mq and effective radius ~ 10 kpc (Buote et 
al.|2004| . 



We start our calculations with the hot gas in spherical 
hydrostatic equilibrium in the potential well generated by 
the total system mass. We use the observed T{r) and n(r) 
( [Buote et aL][2003l [2004} see the dotted lines of Fig. 1) to 
calculate the total gravitational potential under the assump- 
tion of hydrostatic equilibrium. At the virial radius the gas 
fraction is ~ 0.11, a reasonable value for X-ray bright groups 
( [Mathews et al.[[2005[ ). 

The computational rectangular 3D box in all of our 
models extends slightly beyond the group virial radius, iiio4. 
We simulate the z ^ half-space with symmetric boundary 
condition at z = 0, while elsewhere we set prolonged initial 
conditions with only outfiow allowed. Despite the AMR ca- 
pability of FLASH, we decided to use a number of concentric 
fixed grids in cartesian coordinates. This ensures a proper 
resolution of the waves and cavities generated in the cluster 
core by the AGN outfiows. We use a set of 10 grid levels 
(basic blocks of 8 x 8 x 4 points), with the zone linear size 
doubling among adjacent levels. The finest, inner grid has 
a resolution of 488 pc and covers a spherical region of ~ 20 
kpc in radius. In general, grids of every level extend radially 
for about 40 cells. This simple method is the best way found 
to cover large spatial scales (hundreds pc up to Mpc) and at 
the same time integrating the system for several Gyr, using 
moderate computing resources. 



2.2 Outflow generation 

We adopt a purely mechanical AGN feedback in form of 
nonrelativistic, coUimated outflows. The implementation is 
similar to the galaxy cluster outflows of Gil. 

In considering massive slow outflows, we are implicitly 
assuming that the relativistic jet entrains some ICM mass 
(Afact), or in a similar way, that radio jets are highly rela- 
tivistic on pc scale, but rapidly decrease to subrelativistic 



velocities within few kpc from the black hole (Giovannini 



20041 



As in Gil we show results only for models with cylin- 
drical jets, with velocity parallel to the z— axis. We have 
calculated few simulations with conical (or precessing) out- 
flows and verified they have a similar impact on the global 
properties of the flow. In fact, the pressure of the surround- 
ing gas coUimates the outflows within few tens kpc. 

In this paper we will focus mainly on the two most suc- 
cessful models adopted for clusters in Gil: cold and Bondi 
feedback. 

In the cold feedbaclj^an outflow is activated only when 



See also Soker|2006 for a similar type of heating. 
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gas cools to very low temperature within a spherical re- 
gion r < 3 kpc, and drops out from the flow. If AMcooi 
is the mass cooled in the aforementioned region during a 
single timestep, the injected kinetic energy (in the following 
timestep) is given by: 

AEjet = e AMcool (6) 

This energy is given to the hot gas located in a small re- 
gion at the centre of the grid (the 'active jet region'), whose 
size is alway^2 x 4 cells (~1 kpc wide, 2 kpc high), con- 
taining a hot gas mass M^^t- 

\M^ctvl, = ASj,t. (7) 

We will see that the frequency and strength of the feed- 
back events strongly depend on the mechanical efficiency e, 
which has typical values I0~* — 10~^ (see Table 1), almost 
an order lower than galaxy clusters. 

In the Bondi feedback the outflows are triggered by the 
usual prescription: 

Mb = 47r(GA/BH)Vo/cfo, (8) 

where po is the volume-weighted hot gas density calculated 
within ~ 2.5 kpc, while Cso is the mass-weighted sound speed 
in the same region. The assumed SMBH has a mass of 3 x f 0^ 
Mq. The outflow energy is then calculated with eq. (7), using 
Mb instead of AAfcooi. Needless to say, the Bondi radius, 
tb ~ 80 pc, is smaller than our resolution limit, so we refrain 
to attach a strict physical meaning to Mb . In this sense the 
high mechanical efficiencies adopted for Bondi models are 
due to the fact that the accretion should be considerably 
larger (because of higher inner po and lower Cso). 

Besides the two main feedback schemes, we exper- 
imented some variations in the AGN generation, such 
as fixed intermittency, thermal feedback, direct linking 
between outflowing and accreting mass. These special 
models are often not successful and they will be briefly 
explained through Sections 3.4 - 3.6 (see also Table 1 for a 
summary) . 

For most of the presented methods we applied injection 
directly in the domain with Mact, as stated above. In other 
simulations we tested also the injection of mass, momentum 
and energy flux through a 'nozzle' in the grid boundary at 
z = Q. This way, the jet power is expressed as: 

^ (pjotUjot^n)uj^ct = eAf (9) 

with a nozzle area A-n of 2 x 2 cells and pjct ~ O.I the initial 
central gas density. We fix the temperature of the jet to very 
low values, compared to the IGM, in order to keep injected 
thermal energy on negligible levels compared to the kinetic 
flux. 

As noted also in Gil, the results do not greatly depend 
on the method of injection. However, the physical approach 
is a bit different: in the Mact method we are directly mod- 
elling the entrainment, while in the nozzle injection we are 
specifying the properties of the 'subgrid jet'. 



A slightly smaller or bigger jet does not alter the global evolu- 
tion, see Gil. 
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Figure 1. Evolution of CF model (no AGN feedback). Top panel: 
gas cooling rate versus time. Middle and bottom panels: temporal 
evolution of the gas (electron) number density and mass weighted 
temperature profiles, respectively. The profiles are displayed at 15 
different times, as indicated in the lowermost panel. Chandra and 
XMM-Newton observational data of NGC 5044 are represented by 
star and cross points l |Buote et al.|2004[ |, fill ed circles l|Buote et| 
|al.|2003| l; open circles indicate ROSAT data. l |David et al.|1994| l. 

3 RESULTS 

In this Section we report the results of the various simula- 
tions, exploring different types of feedback and parameters. 
We have analysed in detail the long term behaviour of mod- 
els similar to the most successful ones in Gil. Some of the 
properties of unsuccessful - but pedagogical - models are 
covered in Gil. 

We stress again that the main objective of the present 
work is investigating the global properties of the flow, such 
as cooling rates and azimuthally averaged profiles (density, 
temperature, etc.). Cavities, shocks and iron abundances 
will be studied in-depth elsewhere through a dedicated set 
of simulations at higher resolution. 

3.1 Pure cooling flow 

As a fiducial model, we ran a simulation without AGN feed- 
back, i.e. a pure cooling flow (CF). The results are shown in 
Fig. 1. 

The flaws of the classical cooling flow model are clearly 
depicted in this evolution. Both the density and (mass 
weighted) temperature profiles steepen considerably in the 
central region: the radiative losses induce a subsonic inflow 
of the gas, increasing the emissivity in a vicious cycle. This 
trend is relatively more evident and fast in the galaxy group. 
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rather than in the cluster. After just 1 Gyr the coohng rate 
has aheady approached a quasi steady state, with an asymp- 
totic value ~ 25 M0/yr. This is the clearest discrepancy 
with observational data and one of the main failure of the 
standard cooling flow scenario. 

The other two main observables, T{r) and ne(r), are 
also clearly wrong in the inner ~ 15 kpc (see Sec. 1 and 2.1 
for the observational references). Very cold gas, down to 10^ 
K (which is soon dropped in our simulation), concentrate 
in the nucleus and generate non consistent gradients. Our 
feedback models have to reduce this steepening, but preserv- 
ing the cool core appearance at the same time. As said in 
Sec. 1, galaxy groups present a more variable inner temper- 
ature profile respect to clusters. Nevertheless, the heating 



must not destroy the inner structure (see also Brighenti & 



[Mathews, 2002 2003), a quite demanding requirement. 

A typical signature of such a strong cooling flow is also 
the temporal increase of the bolometric X-ray luminosity: 
from ~ 1.9 X 10*^ erg s"^ up to ~ 2.6 x lO''^ erg s~^ at 7 
Gyr. 

It is interesting to investigate the global energetic bud- 
get, because the common (misleading) sense would proba- 
bly associate the irradiated energy with, mainly, the drop 
in internal energy. However, the internal energy within r^ii 
decreases by ~ 7 x 10^* erg (that is, by less of 1%), while 
the potential energy drops by ~ 6 x 10®° erg, considering 
both the hot gas remaining in the grid and the cooled gas 
at the centre of the cluster. The kinetic energy stays always 
around ~ 10^^ erg and is therefore negligible. The conclu- 
sion is that energy is radiated away (-Erad ~ 6 x 10®° erg 
in 7 Gyr) mainly at the expense of the potential energy of 
the IGM. This behaviour was less evident in the cluster, but 
still present. 

In order to test the effect of numerical resolution on the 
CF results, we calculated a CF simulation with a resolution 
three times higher: the results are almost identical to the 
lower resolution run described above. Therefore, with the 
adopted resolution our models are clearly in the convergence 
limit. 

After placing the frame of the requirements for a suc- 
cessful feedback model, in the following we will describe the 
simulations with heating linked to AGN outflows. In the next 
Section, we will begin with the two best Bondi models and 
then we will discuss, from higher to lower efficiencies, com- 
putations with worser results, alternating continuous and 
intermittent type of feedback. 



3.2 Bondi feedback (or entropy-regulated) 

3.2.1 Model Bc5em2, es = 5 x 10"^ 

We begin to analyse the results of the simulations with AGN 
feedback turned on. 

The Bondi accretion theory (Bondi 19521, although 



highly idealized with respect to the complexity expected in 
real accretion systems, has been widely used to estimate ac- 
cretion rates on supermassive black holes, adopting Eq. Q. 

Despite our resolution is far from capturing the effec- 
tive Bondi radius (tens of pc) and despite the fact that the 
radiative cooling is important in the central region, we think 
that Bondi accretion is still a useful, simple feedback pre- 



scription. In fact, the accretion rate is simply proportional to 
K~^^'^, where K = T j (?^^ is the entropy parameter. Thus, 
Bondi regulation is intrinsically a good 'thermostat': when 
the gas is cold and dense, the feedback is stronger, and vice 
versa. 

The first good model that we obtained is Bondi accre- 
tion with efficiency e = 5 x 10~^ (Fig. [2] second column), 
using the 'entrainment' injection in the active region. The 
zone for the weighted quantities is ~ 2.5 kpc. The overall 
feedback is gentle and moderate, exactly as speculated. 

The mean cooling rate is well below 10% of the pure 
cooling flow, with only a transient peak (~ 3.5 Mq yr~^) at 
1.5 Gyr. An important feature for a successful continuous 
feedback, as seen in Gil, is stability in time. The moderate 
power stays around 10** erg s~^, without strong variations, 
which allows to balance the instantaneous cooling rate in an 
efficient way. 

The mass-weighted temperature profiles are in agree- 
ment with observations in the nucleus of the group, with 
a smooth and almost fiat curve. In the outer regions tiny 
acoustic oscillations are present, while between 10-30 kpc 
the continuous injection of energy tends to slightly overheat 
this region, especially at late times. The density profiles also 
show little variation from the initial observed ones, an indi- 
cation that the gas is not accumulating and peaking in the 
core. 

Another relevant consideration is that the Bondi rate 
is always sub-Eddington (between 10~^ — 10~^ Eddington 
rate). Observations of X-ray bright elliptical galaxies (e.g. 
Allen ct al. 2006) also point towards similar values. 

The final injected energy is relatively low, ~ 1.8 x 10®^ 
erg (3.6 x 10®^ erg for the bipolar outflow), compared to the 
total 'available' BH energy 1.8 x 10®^ erg; roughly estimated 
as E'er ~ 0.1 Mbh (? , with Mbh = 10^. This wiU be a typ- 
ical feature of galaxy groups, in contrast to the high input 
energies of massive clusters, whose models can sometimes 
exceed the above quantity. 

Furthermore, we tested several initial conditions and 
we concluded that the dynamics and results of our feed- 
back models are, overall, not distorted by relatively differ- 
ent initial temperature and gravity proflles. For example, in 
Fig. |2] (fourth column), we started the computation with an 
isothermal profile (~ 1.35 x lO'^ K), retrieving the density 
structure from a theoretical NFW plus deVaucouleurs po- 
tential. After the first Gyr of pure CF, the Bondi feedback 
(eb = 5 X 10~^) is again very effective in quenching the 
cooling fiow and at the same time preserving the cool core 
appearance. 

A possible riddle for all continuous models remains the 
absence of frequent jet-infiated spherical cavities. The con- 
tinuous AGN activity (mean v ~ 10000 km s~^) carves a 
narrow tunnel of about 30 kpc in length, although its den- 
sity contrast with the environment is large only for z < 10 
kpc (and not particularly evident in the surface brightness 
maps; see Sec. 4). 



3.2.2 Model Bi5em2, eb = 5 x 10"^ 

In a self-regulated (subgrid) mechanism we have to face sub- 
stantially two problems. The first one is how much power to 
link to the feedback, while the second one is the right timing, 
i.e. when to turn on and off the heating machine. 
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Figure 2. Evolution of Bondi (continuous) feedback models with different efficiencies (tfie latter being the isothermal start). Panels 
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kinetic energy (in the half-space z > 0), outflow power and velocity, associated cold efficiency (see text). The label 'xlO' means that the 
units in the ordinate axis have been multiplied by 10. 
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In the following simulation (Figure [s] third column) we 
link the activation of the AGN heating to the gas cooling, 
but use the standard Bondi rate for the gas accretion onto 
the black hole. In other words, an outflow is generated when 
AMcoid 7^ 0, with an energy linked to Mb in the usual way. 

This feedback scheme tries to reproduce a self-regulated 
intermittency of the outflows (ideal Bondi models can only 
be continuous), with cycles of inactivity corresponding to 
moments of high hot gas entropy. In fact, when the gas be- 
gins to cool, it will usually start to flow towards the central 
black hole. 

This way, the temperature profiles are a bit more consis- 
tent with the observational data, especially at intermediate 
radii. Although the jets are intermittent, their frequency is 
high, a condition that grants the quenching of cooling rate 
{Mcooi ~ 0.5 — 2 Mq yr~^). The density profiles do not de- 
part from the cool core status, without a strong steepening 
(central values ~ 0.1 cm~^). 

Bondi accretion rates are almost identical to Bc5em2, 
with the total energy injected in 7 Gyr being a bit less: 
~ 1.5 X 10®^ erg. The average jet power is also a factor of 
~ 1.3 less, ~ 5.5 X 10"*^ erg s'^ 

This model clearly solves the problem of jet-inflated 
cavities in the IGM. The onflows duty cycl^is roughly 80-85 
per cent (in the whole computation). The pattern is not very 
regular, in fact the cycle is more frenetic between 2 and 5 
Gyr, when the central gas has been heated efficiently. Every 
time the jet (with velocities usually of 8000-10000 km s~^) 
turns on again, a bubble of small size (10 — 15 kpc) is seen 
in the density map. Note that after a period of quiescence 
the subsequent injection produces a peak in the power, pro- 
moting the above mentioned cavity inflation through shock 
heating. The side effect consists in more evident perturba- 
tions in some T profiles (e.g. 3 and 7 Gyr), however still 
compatible with typical galaxy group observations. 

Finally, during runtime we calculated, at every 
timestep, the efficiency (ec) that would have the cold feed- 
back mechanism (bottom panel). The result is striking: on 
average, the two schemes are almost coincident with ec 
around 10~^, except for few short timesteps, in which the 
cold feedback requires a 100 per cent efficiency (because the 
cooled mass is almost zero). 

3.2.3 Model Bcleml, eb = 10"^ 

We begin now to describe models from higher to lower 
efficiencies, alternating again continuous and intermittent 
Bondi methods. 

The following simulation is very computationally expen- 
sive, because the (continuous) outflow velocity rises often 
over 20000 km s~^ It is shown in Fi gure 2, third column. 
The cooling rate is quenched to fraction of M© yr~^ after 
just a few Gyr. The total injected energy is 3 x 10^^ erg 
(about 1.5 times greater the intermittent counterpart, next 
model Bileml). The consequence of this high energy deposi- 
tion is a progressive steepening in the temperature gradient, 

^ The common duty cycle retrieved by observations can be quite 
different from our numerical estimates, because the definition of 
an 'active' AGN highly depends on some luminosity threshold (in 
the X-ray, radio, etc.) and on the properties of the sample. 



with a moderate spike in the inner 4 kpc, where the jet is 
always active. During some timesteps the outflow reaches 
5 X 10** erg s~^, a quite strong power for Bondi accretion. 

The compared efficiency analysis (bottom panel) shows 
that a cold feedback triggering could not even be possible in 
20-30 per cent of the evolution, because it would break the 
unity threshold. In the other periods it would be associated 
with a value around 2.5 x 10~^. 

3.2.4 Model Bileml, eb = 10~^ 

The following (Fig. |3] fourth column) is another computa- 
tionally expensive model, which touches the limit of accept- 
ability in the parameter space, for Bondi models. The timing 
is again associated with the cold gas. The cooling rate is very 
similar to the precedent continuous model (Bcleml), with 
the reasonable range of 0.5-1.5 M© yr~^ range, declining at 
late times. The density profiles are optimal, not altered from 
the initial data. Temperatures start to feel the strong out- 
bursts, which have sometime power in excess of 10*^ erg s~^ 
(with velocity larger than 2 x 10* km s~^). In fact at later 
times the central gas is heated up to 10^ K, during some 
event. 

The associated cold efficiency has a common trend 
around 2.5 x 10~^, but with many oscillations after 3 Gyr, 
reaching easily unity. Notice that these moments last very 
few timesteps, after which reset to lower values (in con- 
trast to a first impression of the tight plot). 

3.2.5 Model Bclem2 and Bilem2, es = 10"^ 

When the efficiency is lowered to e = 10~^ the heating pro- 
vided by the feedback is considerable (~ 8 x 10^" erg). This 
decreases the cooling rate down to at least 30-40 per cent 
of that of the pure CF (Figure [2j first column). The profiles 
are close to the observed data, thanks to an always contin- 
uous outflow with velocity 4 — 5 x 10^ km s~^ and power 
4 - 5 X 10*^ erg s'^ 

The intermittent simulation is very similar to the con- 
tinuous evolution (Fig.js] second column), except that, ini- 
tially, the system is permitted to cool. Without this early 
opposition, the outflows - with same power - can not halt 
the cooling so efflciently like in the previous case. The feed- 
back finds an equilibrium only at 4 Gyr, while at later times 
it will slowly increase, on the contrary of a reasonable self- 
regulated model. 

The associated efficiency of a AMcow model is near 
10~*. In fact, we will see that the linked cold feedback model 
will produce an almost identical outcome, i.e. weak heating 
(Fig. |4j second column). 

3.2.6 Model Bilem3 (or BclemS), es = 10"^ 

As a template for a typical weak model, a simulation with 
Bondi efficiency of 10~^ produce a cooling rate, which is 
identical to a pure CF run: ~ 25 - 26 Mq yr"^ (Fig.jsj first 
column) . 

The only clue for the presence of AGN heating is that 
the profiles do not suffer a heavy decline with time, even if 
they posses a remarkable steep gradient (T down to 4 x 10^ 
K). 
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The outflow is always continuous (the model could be 
also named BclemS), with a very low velocity around 1000 
km s~^ and a power often under 10^'' erg s~^. There is no 
possibility that this model will become intermittent at any 
time, due to the constant central cooling. 



3.3 Cold Feedback (AMcooi regulated) 

3.3.1 Model C5em5, ec ^ 5 X IQ-^ 

The results presented in the previous Section for the Bondi 
accretion models suggest that, compared to galaxy clusters, 
groups require lower efficiencies of at least a factor of 10. 
Therefore, we start illustrating the cold feedback mechanism 
for the lowest adopted ec and then moving to more efficient 
outflows. 

Here we still use the entrainment injection method, with 
width and length of the active region 1 and 2 kpc (2 and 4 
grid points), respectively. In Gil, in fact, we concluded that 
the size of the jet and the type of injection do not change 
drastically the evolution of the feedback. 

As shown in the first top-left panel of Figure |4j the 
cooling rate settles on a steady value around ~ 10 M0 yr~^. 
This rate is obviously unacceptable, although reduced to 
the ~ 50% of the pure CF model value. With such a low 
efficiency, the azimuthally averaged, mass weighted temper- 
ature and density profiles are only slightly modified by the 
heating and do present the excessive accumulation of gas 
near the center, typical of the CF model. In the center the 
temperature is ~ 5.6 x lO'^ K and the numerical densities 
~ 0.18 cm"^. 

The outflows, typically have a power of 4 — 5 x 10^"^ erg 
s~^ and a velocity of 3000 — 4000 km s~^ . The total injected 
mechanical energy is in fact quite low, 6 x lO^'' erg, and the 
outflows are almost continuous, because they are never able 
to stop the cooling flow. Given the rare period of inactivity, 
the 'duty cycle' is ~ 95 per cent. 

At the end of the simulation 6 x Mq have cooled 
and dropped out of the hot phase. The vast majority of gas 
cools at the very center. If all the cooled gas were accreted 
onto the central black hole, as we have simplistically as- 
sumed in this scheme, the final BH mass would result about 
an order of magnitude greater than the expected one. In 
principle, we could avoid the problem of the excessive black 
hole mass by assuming that only a fraction of the cooled gas 
actually accretes on it with a higher efficiency (see Sec. 3.6, 
InOut models). 

3.3.2 Model Clem4, = 10"* 

As anticipated, the next step is to increase the efficiency, in 
order to find a more suitable heating power. 

In the following model, where the efficiency increased 
by a factor of five, ec = 10"'' we obtain in fact a jet power of 
« 10** erg and velocity oscillating around 5000 km s~'^ 
(see Fig. [4] second column). However, the cooling rate is only 
slightly reduced to 6 — 7 Mq yr~^ and radial profiles tend 
to be similar to the ones of previous model. The difference 
is that now the curves decline at 6 kpc, instead of 10 kpc. 
Evidently, the outflows deposit their energy in more inner 
zones. 

We conclude that, even if the total injected energy 



reaches almost 10 erg, at this stage the cooling flow is not 
quenched in a successful manner. This is also underlined 
again by the almost continuous presence of the outflows. 



3.3.3 Model C5em4, 



5 X 10" 



At 5 X 10~* we finally recover a successful heating evolution. 
As in previous computations, at early time (first Gyr) the ra- 
diative cooling prevails over the outfiow mechanical heating 
and the cooling rate reaches a peak of 4 Mq yT~^ . However 
the long term evolution presents different intriguing results. 
After 2 Gyr the system becomes much more turbulent and 
chaotic, inducing a reasonable amount of mixing, that per- 
mits the gas circulation and reheating in a very efficient way, 
also in the 5 kpc nucleus. 

As shown in the third column of Fig. |4j the temper- 
ature in that region is almost constant, with values close 
to that observed for NG C 5044 (r; 8 x 10'^ K; |Buote et alT 
2003 David et al. [2009 1. This very important feature is also 
described by the density curves: a gradual flattening of the 
central profile is seen through the entire evolution. It is strik- 
ing that the 7 Gyr (black) profile is practically superposed 
to the initial fit. At that time (redshift — 0) the cooling rate 
is ~ 2 Mq yr^^, about 7% of the pure CF run, a significant 
suppression which brings the value of Mcooi in reasonable 
agreement with the observations ( David et al.||2009 Buote 
et al.|2003||Tamura et aL]|2003[ ). 

The shock waves generated by the outflows, with power 
usually between ~ 2 — 5x 10** erg s~^, are not strong enough 
to significantly alter the global positive temperature gra- 
dient, although small amplitude ripples are present. These 
waves correspond to weak shocks associated with the jet 
propagation and are visible up to a distance of ~ 300 kpc. 
They are not so evident as in galaxy clusters because the 
energy involved here are two or three orders of magnitude 
lower. In fact the total jet energy is 1.8 x 10^^ erg. The gas 
is typically ejected with Vjct ~ 10* km s~^. 

Furthermore, an important difference between clusters 
and groups is the duration of each event. In cluster mod- 
els, a duty cycle of 0.1 was quite usual, because the very 
short super-Eddington bursts could stop instantly the cool- 
ing flow, with the consequence of perturbed proflles, near the 
ignition time. Here, our estimated duty cycle is roughly 0.85. 
Therefore the situation is reversed: successful cold feedback 
models in groups require an almost continuous heating, with 
short pauses of tens Myr. 

Notice that the Eddington luminosity, Z/Edd ~ 1.5 x 
1O^**(Mbh/M0) ~ lO*'' erg (for a 10^ Mq black hole), is far 
above the regime of the current model. In fact the accretion 
rate oscillates between 10"'' — 10~^ the Eddington rate. As 
seen in Section 3.4, a I/Edd jet can easily erase the cool core 
structure of the group. 

In summary, the analysed behaviour of C5em4 seems to 
be similar to that the quiet Bondi feedback. In the bottom 
panel we have calculated, at every timestep, the associated 
Bondi accretion rate and the required eb to reproduce the 
same instant mechanical energy. As a striking result, a Bondi 
efficiency of ~ 5 x 10"'^, especially after 2 Gyr, is retrieved 
from this analysis. Only in a few events, the cold mecha- 
nism detaches from the regular regime of Bondi, because of 
its intrinsic impulsivity. Nevertheless the similitude is quite 
evident and indeed, if we compare the two outcomes of both 
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Figure 4. Evolution of all cold accretion models with increasing efficiencies, from left to right. The description of the plots are same as 
for Fig. 2. 
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models (especially intermittent Bondi, Fig.|8] third column), 
they appear deeply connected also in a quantitative manner. 

3.3.4 Model Clem3, Ec = 10"^ 

Increasing the efficiency by a factor of 2 generates another 
positive model, with some features approaching the border- 
line of a violent disruptive heating. 

In fact, with a jet power up to 10*^ erg s~^, the cool- 
ing flow is perfectly stifled, asymptotically decreasing to a 
fraction of M© yr~^ after a few Gyr (Fig.H fourth column). 

The beginning of the violent regime is quite evident in 
the averaged proflles, which now tend to oscillate in a large 
~ 50 kpc region. Overall, the cool core is conserved, but 
the temperature indicates some overheating in the nucleus 
at later times. 

With our adopted active region, in a few outbursts, the 
velocity is greater than 10^ km s~^ line, while the outflow 
mass rate is several tens Mq yr~^, as in all previous AMc 
models. This is a clear indication that at higher efficiencies 
we will approach a relativistic regime and a catastrophic 
heating (we will not show here results for ec > 10""^ but 
see fig. 5 of Gil for an analogous model for clusters). The 
equivalent Bondi efficiency for model ClemS is sometimes 
0.5 or more, difficult to justify. 

The total injected energy of ClemS is similar to the 
previous successful model, a sign of a good self-regulation in 
the global feedback process. This energy can be compared 
with the total energy radiated away, « 1.38 X 10^" erg (in 
the simulated half-space z > Q). Interestingly, the present 
heating provides more energy than that lost by radiation 
(^-rad ~ 1.4 X 10^'^ erg). Notice that, in every described 
model, the total outflow energy (for the full-space system) 
is about an order of magnitude lower than the available BH 
energy (~ 1.8 x 10®^ erg). The evolution of the energetics is 
very similar to the description presented by Gil. The core 
kinetic energy increases after every intense AGN outburst, 
but it is soon dissipated and transformed into potential en- 
ergy through the expansion of the IGM. 

3.4 Intermittent feedback 

3.4.1 Model IntS 10m: 10'^ PEdd 

In this series of models we want to decouple the simulation 
from the self-regulation engine, in order to control some im- 
portant AGN outbursts parameters. For example, it is inter- 
esting to fix the frequency of the outfiows. In Gil we found 
that intermittent outfiows, with jet duration of 5 Myr and 
activated every 10 Myr, are very effective in inhibiting gas 
cooling preserving the correct temperature and density pro- 
files in galaxy clusters. This AGN activity cycle agrees with 
observations of acoustic 'ripples' in Perseus cluster (e.g., 
Fabian et aL]|2006| ). 

In the models described in this Section we do not insert 
an ad hoc velocity, but we retrieve its value from setting a 
constant power, here taken equal to 10~''PEdd. Thus Pjet ~ 
4 X lO''^ erg s~^. We choose the 'nozzle' injection method, 
therefore the velocity is constant: 7000 km s~^. Note that 
the efficiency is irrelevant in the computation. 

As a result (Fig. |5] first column), the cooling rate is 
less than a few Mq yr~ , for several Gyr, and shows a slow 



increase after 4.5 Gyr, up to 4 Mq yr~^. The azimuthally 
averaged density profiles are almost superposed to the ob- 
servations, with a slight decrease in time, becoming flatter. 

Thus, like in Bondi models, the radiative cooling is sig- 
nificantly reduced by the simulated heating process. Using 
the nozzle injection (essentially an outfiow with negligible 
height) produce a spike in the very few kpc, while at radii 
greater than 10 kpc the fluctuations from initial T{r) are 
almost zero. As expected, shock heating is more vigorous in 
the nucleus when the outflow is injected through the bound- 
ary, while more effective at intermediate radii, when the en- 
trainment method is adopted. 

In this type of intermittent models it is not granted 
that the feedback produces evident cavities. In fact with a 
cycle of 5-10 Myr the head of the jet cocoon almost touches 
the backfiow tail of the last bubble. Thus, we expect that 
with a duty cycle greater than 50 per cent the jet becomes 
practically continuous. 



3.4.2 Model Int510h/l: 10"^ and lO'^^PEdd 

Keeping the same cycle scheme (5-10 Myr), but increasing 
the power by a factor of 10 produces a typical 'exploding' 
simulation (not shown). The temperature profile gets a very 
steep negative gradient, as soon as the outfiow ignites, de- 
stroying the entire cool core and generating a 'plateau' in 
n(r). The cooling rates are obviously zero, but with unac- 
ceptable observables. 

On the other hand, with a much lower power (4 x lO*'^ 
erg s~^) Mcooi approaches 12 Mq yr~^ after just 1 Gyr. Gen- 
erated bubbles are quite stable and, being slowly infiated, 
show cold rims in the early phases. Nevertheless, the den- 
sity and temperature profiles resemble those of a classical 
cooling ffow, with less peaked gradients. 

We can conclude that, even with fixed intermittency, 
an acceptable AGN feedback model should not approach an 
Eddington regime (like proposed in some analytical works, 
e.g [King 2009 ), at least in galaxy groups, where a strong 
burst can easily wreck the delicate thermodynamic structure 
of the core. 



3.4.3 Varying fixed cycle: 1-10, 2.5-10, 3.3-10, 7.5-10 

Keeping a fixed jet power, as in our acceptable model 
Int510m, but this time varying the duration of the AGN 
activity, leads to models that must be rejected. 

Indeed, when every outfiow event has a duration of 1, 
2.5 or 3.3 Myr the intermittent feedback becomes quite inef- 
fective in halting the inff owing cold gas: the central gas den- 
sity slowly grows with time until radiative cooling prevails 
over heating, causing the cooling rate to surpass the thresh- 
old of acceptability. For example after just 1 Gyr, with 3.3 
Myr duration, the cooling rate is over 7 Mq yr"'^ and the 
temperature curve turns very similar to those of model GF. 

On the contrary, when the jet duration exceeds 5 Myr 
the inffated cavities disappear, because the different out- 
bursts melt rapidly in a ten kpc zone, producing an almost 
continuous injection. The consequence is a steep T gradient 
at the center, which can halt the cooling flow, but without 
self-regulation it easily overheats the system. 
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3.5 Thermal feedback 

3.5.1 Adding thermal energy: 20% up to 99% 

It is worth noting that one difficulty of a pure kinetic outflow 
(often caUed 'momentum-driven') is that it can easily carve 
a tunnel in the surrounding IGM if the ram pressure of the 
jet {pjVj) is much larger than the thermal pressure (pc^/'y) 
of the gas. Thus a feedback with high kinetic energy will 
advance undisturbed in the group core, depositing its en- 
ergy at intermediate radii (15 — 30 kpc). At the base, the 
cold gas accumulates in a torus-like manner. This is espe- 
cially evident at early stages of a momentum-driven outflow. 
The successful models, we found, are able to quench cool- 
ing flow, because at some point the generated turbulence 
induces strong mixing in the gas, also in the very center. 
Therefore the circulation inside the jet path makes shock 
heating very effective. In our models the turbulence is sus- 
tained only by the AGN feedback. It is likely that this pro- 
cess dominates in the very central region of the system, al- 
though turbulent and bulk motions caused by cosmological 
accretion, usually very subsonic, may also contribute (e.g., 
Heinz et al.|2006 |. 

Nevertheless, the high density gas in the torus cools 
efficiently and triggers the AGN feedback in an almost con- 
tinuous way. It is physically reasonable that a part of the 
jet power is converted in thermal energy, especially during 
the entrainment process. Based on these considerations, we 
simulated some models with a fraction of the total injected 
energy in the form of thermal energy, through the nozzle. 
We kept the Bondi prescription for the feedback method. 

As pedagogical purpose, we first set the fraction of Eth 
as 99 per cent the total energy. This setup is similar to the 



model of Cattaneo & Teyssier (20071 with loading factor of 



100 and Bondi efficiency 0.1. As expected, the first powerful 
outflow (over 10** erg s~^) destroys the entire thermody- 
namic structure of the cool core: the produced shock (with 
temperatures above lO^'' K) is almost spherical. This is an- 
other indication that the group necessitates a much more 
delicate heating compared to a cluster. Lowering the effi- 
ciency to 10""^ still generates a similar catastrophic heating. 

On the other hand, keeping the above efficiency, but re- 
ducing the fraction of Eth to 20% or 33% results in models 
overall similar to purely kinetic ones. One major difference 
are small puffs at the base of the jet, that enable more tur- 
bulence. The problem is that, in the inner 8 — 10 kpc, tem- 
perature profiles show a spike, due to injected Eth. Another 
interesting effect is the enlargement of the tunnel carved 
by the outflow, because of the increased internal pressure, 
compared to the momentum-driven one. 

The best model seems to be Eth50 (not shown), with 
50% thermal energy and efficiency 5 x 10~^. The feeble jet 
is almost always continuous, with velocity 3000 km s~^, but 
the thermal energy injection generates a break-up in the jet 
structure, which fragments in small buoyant bubbles of size 
5 — 8 kpc; a very interesting feature that is seen in NGC 
5044 dPavid et al.|2009| [2010 K 

We conclude that a good fraction of thermal energy 
associated with the AGN feedback can help the deposition 
of heating at very small radii, reducing the cold torus at 
the base of the jet and fragmenting the outfiow in small size 
bubbles. At the same time, however, it is difficult to obtain 



a flat or positive (mass-weighted) T profile. We will expand 
this analysis in a future work, with very high resolution. 



3.6 In & out feedback 

3.6.1 Model 101-9 

One problem of self-regulated feedback models is how much 
of the gas contributes to the growth of the black hole mass. 
Needless to say, with a resolution of ~ 500 pc we cannot 
consistently estimate the accretion rate onto BH. As noted 
in some works (e.g., Cattaneo fc Teyssier||2007 Ostriker et 
ar||2010 Gil) part of inflowing gas onto a black hole could 
be thrown back by the generated outflow, perhaps through 
entrainment. Therefore only a fraction of the accreting mass, 
AMacc, might actually be captured by the black hole. From 
another perspective, we can see this as a reduced total effi- 
ciency. 

That being stated, we computed some models in which 
we attempt to track the part (/in) of infalling gas that really 
increases the BH mass. The residual fraction (/out = 1 — /in) 
is considered the active mass of the outflow: 



-/out Ma 



c^j^ct = /in eMa.ccC . 



(10) 



Note that in these models it is better to use the nozzle injec- 
tion mode, in order to totally control the outflowing mass, 
without altering the internal domain. 

An interesting feature of this method is that the out- 
flow velocity is practically constant, dependent on /in//out 
(and e) . An observational constrain of the outflow velocity 
may thus give hints about this ratio between infalling and 
outflowing matter. 

Setting /in = 0.1 (/out = 0.9) and efficiency 10"^, the 
outflow, coming out of the nozzle, has a velocity around 7000 
km s~^. With cold accretion, the power of the jet presents 
strong oscillations between 10*^ and 10*^ erg s"'^ while the 
cooling rate approaches in a few Gyr the pure CF simu- 
lation. Even with Bondi accretion (same efficiency; Fig. |5j 
second column), which, in theory, is more regular, the (high) 
outflow power can not find a stable balance in time, result- 
ing in Mcooi ~ 18 Mq yr~^ at 7 Gyr. While some accretion 
events reach 4-5 M© yr~^, the produced strong feedback 
does not last long enough. In fact the total injected energy 
is very low: 10^^ erg. It is interesting that the radial profiles 
of this model (101-9) show a cool core, with the exception of 
a central negative temperature gradient. Therefore, overall 
acceptable mean T and n profiles do not necessary imply 
that the system is not cooling. 

Trying to increase the efficiency is a serious problem 
for the InOut model, because the velocity will constantly 
overtake 50000 km s~^. We could lower /in in order to in- 
crease the jet mass, but with just 1 per cent of accreting 
mass (probably unrealistic) the velocity will still have very 
large values. 



3.6.2 Model lOeniO and lOenSO 

The only solution, we have found, is to return back to the 
'entrainment hypothesis' and multiply /out by a factor r] 
(like a mass loading factor). Following this assumption, with 
?7 = 40 and /out = 0.3, we could enhance the Bondi efficiency 
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(e = 5 X 10""^) to obtain an acceptable cooling rate, below 
5 Mq yr~^ (Fig. [5] third column), after the small bump at 
1.5 Gyr. Now the velocity has a value of ~ 6000 km s~^. 
With = 80 the average cooling rate is further reduced, 
asymptotically decreasing to 1 yr~^ (fourth column). 
Both simulations generate some of the best radial profiles, 
which keep the cool core appearance for the entire evolu- 
tion time, almost 'glued' to the initial (observed) condition. 
This is an indication that the feedback has returned stable 
(around Pj ~ 1 — 3 x lO**^ erg s~^), gentle and continuous, 
very similar to successful Bondi models (Sec. 3.2). Notice 
that the accretion rate is again sub-Eddington. 

Overall, a pure InOut model exhibits a spasmodic be- 
haviour, because of the intrinsic linking between accreting 
and outflowing matter, which can be smoothed out only by 
introducing a mass loading factor. However, this parame- 
ter is nowadays unknown and therefore its value is chosen 
ad hoc, in order to retrieve reasonable jet velocities. It is 
therefore simpler and cleaner to adopt the usual entrain- 
ment feedback of the main models (Sec. 3.2 - 3.3), seeing 
that the results are analogous. 



4 DYNAMICS AND OBSERVABLES 

In this Section we will show the 'real face' of the best simu- 
lated models, at some relevant times. We will analyse phys- 
ical quantities, such as gas density, and other astronomi- 
cal observables, i.e. integrated emission-weighted maps. No- 
tice that the entire structure of our computation favours the 
study of a very long term evolution. In order to fully com- 
prehend the dynamics of single burst events, we should carry 
on a short term evolution with much higher resolution. We 
will confront this setup in a future dedicated paper. 

In Fig.[6]the first column exhibits two-dimensional den- 
sity cuts trough the x-z mid-plane (panels a-d) with the 
velocity field superposed. The first row is associated with 
Bc5em2 model (at 5 Gyr), one of the successful models in 
our set of simulations. This snapshot depicts one of the three 
typical stages of the AGN-regulated simulation in a galaxy 
group. The gentle and continuous heating of the outflow 
(~ lO** erg s~^) carves a narrow tunnel in the intragroup 
medium with electron density around 2 x 10~* cm~^. The 
diameter of the channel is usually 4-5 kpc, coUimated in the 
inner region by the thermal pressure of the surrounding gas. 
At 14 kpc, where the jet ram pressure becomes comparable 
to the thermal pressure of the ambient gas (Pth ~ 2.5 x 10~^^ 
erg cm~^), the outflow structure becomes perturbed and 
more turbulenl[^ In fact, near the tunnel, instabilities are 
visible, a signature of further IGM entrainment and turbu- 
lent mixing. The latter has a relevant role in the feedback 
process, especially at later times, when the AGN has been 
active for several Gyr (Fig. [2j cooling rate decreases with 
time). In fact, even if the turbulent mixing does not directly 
heat the gas, it greatly favours the deposition of energy at 
the base of the jet, where the cooling tends to dominate. 

The fact that the global cool core appearance is not 
greatly modifled, is another key result. Although the contin- 
uous outflow generates a tunnel, the overall initial structure 



^ At 24 kpc the outflow is almost stopped by the IGM, generating 
a rise in temperature, ~ 1.55 X 10^ K (see also Fig.[7]l. 



is conserved. On the other hand, increasing the efficiency of 
Bondi models leads to a wider V-shape channel, which be- 
gins to dominate the core of the group. This is not observed 
and in fact those models are rejected also for the central neg- 
ative T gradient. As seen in the second column (panel e), 
the (emission-weighted) temperature of Bc5em2 has instead 
a positive gradient, almost flat in the 5 kpc nucleus. 

The continuous presence of the channel (interrupted 
sometimes by the fragmentation) is a feature not commonly 
seen in real groups, at least in the local univers^ However, 
in the third column (panel i) the X-ray surface brightness 
map does not show an evident channel, with only two very 
faint features (at 12 — 16 kpc), similar to a pair of 'arms'. Go- 
ing into quantitative details we performed one dimensional 
(ID) cuts through z of SBx at different levels (from 4 up 
to 20 kpc; see Fig. [t] non-black lines). It is clear that the 
central depressions are not deep (15 — 20%, usually with a 
difference of 5 x 10~^ erg s~^ cm~^), plus the width is par- 
ticularly slim (few kpc), making the detection of the channel 
extremely difficull[^ Over 20 kpc the slim tunnel is lost in 
the integration trough line of sight. 

The last column is associated with the tracer of the iron 
abundance (again emission- weighted), injected by SNIa and 
stellar winds of the cD galaxy (see Gil for details). This 
type of advected quantity is deeply linked to the history of 
the AGN heating. In panel o is indeed possible to recognize 
the outfiow pattern, which is asymmetrical and can drag the 
metals produced deep inside the central elliptical gala:xy up 
to 15-20 kpc, with values around 0.3 — 0.4 Z©. This char- 
acteristic is commonly observed in the core of AGN heated 
systems (e.g 



|Ponman||2009 



tion) . 



Kirkpatrick et al. 



2009 2011 Rasmussen & 



David et al.||2010| Doria et al., in prepara 



The second row of Fig. |6] shows another good model, 
one based on the cold feedback. Clem3 (at 0.5 Gyr) shows 
the other face of the AGN-regulated evolution: an AGN 
burst with power ~ 1.2 x lO'*^ erg s~^ generates a big cavity 
in the IGM density, with major/minor axis of 16/10 kpc. 
The injected energy is almost 10®" erg, with an outfiow ve- 
locity over 10** km s~^. The buoyant bubble has a high den- 
sity contrast with the environment, ~ 100. 

In contrast to the violent cavities generated in our pre- 
vious galaxy cluster simulations (Gil), the cold accretion 
mechanism is able to inflate low density cavities without 
heavy shocks. The outflow power is in the group simula- 
tions several orders of magnitude lower and the injection 
per timestep is usually smoother. In fact (see panel f and 
Fig.]?! black line of bottom panel), ~ 1.45 x 10^ K inside 
the shell which is only slightly larger than the surrounding 
gas value. This is also suggested by the homogeneity of the 
map above 12 kpc (reddish colours). 

A striking feature, usually hard to detect, is the bright 
rim, that surrounds the cavity (panel 1). The rim, especially 
the low-z region, is formed by low-entropy gas originally 
in the center. This feature is often seen in deep X-ray ob- 



servations (e.g. Salome et al. 2006 and references therein). 



^ The evolution of groups is still far to be certain, lacking a com- 
plete sample at intermediate and high redshift. 
* In addition, Poisson noise, background, resolution, and re- 
sponse of the X-ray detector (not present in our mock maps) 
could easily alter or obscure this very faint and narrow feature. 
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Figure 6. Maps of four best models (kpc unit). Columns, from left to right: cuts through x-z mid-plane of electron number density 
(cm""^, logarithmic) with velocity field superimposed; emission-weighted temperature maps (K, logarithmic); X-ray surface brightness 
maps (erg s~^ cm~^); emission-weighted iron abundance maps (Zq unit). Rows, from top to bottom, are associated with different models: 
Bc5em2 (5 Gyr), ClmeS (0.5 Gyr), Bi5em2 (7 Gyr), Bileml (6.5 Gyr). The colour scale is given by each bar at the top. 



suggesting that the rims are colder than the average ambi- 
ent medium. The ID cut through z = 20 kpc (Fig. [tJ clearly 
confirms the presence of this kind of rim, a drop of temper- 
ature (about 40%) coincident with the high X-ray emission 
(at X = ±6 kpc). 

In panel p we show the (emission-weighted) iron abun- 
dance map in the x — z plane. The iron-rich core, few kpc 
in size, is clearly visible. At z ~ 28 kpc there is a region 
of Fe-rich gas, lifted by the outburst. The dense cavity rims 
also have a slightly larger abundance than the ambient gas. 



also revealing that the origin of a part of the rim material 
is connected to the nucleus of the group. 

The SBx inside the cavity (panel 1) exhibits a depres- 
sion of ~ 50 per cent with respect to the rims. Notice how 
this feature (Fig. [t] black line) would be easily detectable in 
the X-ray, compared to the vanishing faint channel, due to 
its deep depression and large width. Moreover, the nucleus 
(5-8 kpc) dominates the emission, while the upper part of 
the bubble vanishes rapidly in the background (> 30 kpc). 
The future of this bubble is to buoy outwards and being 
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destroyed, after few tens Myr, by the backfiow and instabil- 
ities. 

The third row presents the snapshots of model Bi5em2 
at 7 Gyr. At that stage of evolution the IGM has accu- 
mulated more turbulence, after several AGN events, pro- 
moting the diffusion of iron in a radius of 20 kpc from the 
center (panel q). However, this particular moment is again 
dominated at the bottom-center by another AGN outburst 
(2x 10^'' erg s~^). Therefore, the young bubble (panel c) still 
buoys undisturbed in the IGM, generating a slight asymme- 
try in the iron distribution along jet axis, over the more 
uniform background. 

At this time we caught indeed the early injection phase 
of the outflow. A qucisi cylindrical cocoon envelops the core 
of the jet, which has now a velocity around 3000 km s~^ . The 
injected initial velocity was over 10'' km s-\ meaning that 
the jet commonly decelerates in a rapid way after few kpc. At 
z = 22 kpc the emission-weighted temperature reaches the 
maximum value of 1.62 x lO'^ K (panel g), clearly indicating 
the effect of shock heating. Around the contact discontinuity 
the gas is weakly shocked, with temperatures above 1.25 x 
lO'^ K. The average Mach number is around 1.2, a typical 



Gitti et 



value found by observations (Blanton et al. 2009 
al.|2010 1. The arc is also visible in the SBx map (panel m) 



while the decrement associated with the cavity is ~ 35 — 40 
per cent. 

Notice that the maps of this model are very similar to 
those of C5em4, another successful model, which has anal- 
ogous bubbles production. This behaviour confirms further 
the strict relationship between (intermittent) Bondi and cold 
models, as previously noted by the compared efficiency anal- 
ysis. 

The last row of maps displays another snapshot, typ- 
ical for intermittent (self-regulated) feedback models, like 
Bileml. The big outburst and the generated bubble have 
been vanished, leaving the system in the aforementioned tur- 
bulent phase, in which mixing distributes the low residual 
heating to the central IGM. Soon after this period, the gas 
will cool and start the inflow. The velocity field is relatively 
chaotic (panel d), with maximum velocity of 640 km s~^, a 
bit lower than local speed of sound, evidently decreasing in 
the turbulent zones. The mixing of gas at this epoch restores 
the spherical symmetry of the cool core, with average tem- 
perature of ~ 9 X 10® K (panel h). This is underlined also by 
the symmetric surface brightness map (panel n) and by the 
substantial diffusion of Fe in the intermediate zone, 10 — 30 
kpc (panel r). Note that the abundance radial profile at the 
very center is, however, more peaked during this quiescent 
cooling flow phase (~ 0.8 Z0), compared to the AGN out- 
burst period. Thus, our AGN models might also explain the 
dichotomy of cool core and non cool core groups, having the 
former higher central abundances ( Johnson et al.|[2011 l. 

An interesting feature, captured at this stage (often 
present after an AGN outburst), is the sequence of very 
small 'bubbles' at r > 15 kpc, with typical diameter of a few 
kpc. The origin of this tiny 'bubbles' is the fragmentation 
of the jet due to turbulent and chaotic motions, previously 
discussed. In this map they are in a very late phase, almost 
disappearing in the ambient medium. In this case they 
would be probably not detected in the X-ray (just a 10 per 
cent jump in surface brightness). To capture their signature, 
an observation should be done at an earlier evolutionary 
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Figure 7. X-ray surface briglitness and emission-weighted tem- 
perature ID cuts through z (see legend). Black thick lines corre- 
spond to the model Bileml (panels n and h of Fig. |6|, while the 
others belong to model Bc5em2 (panels i and e). 



moment, after a big event, when the jet is rapidly collapsing. 
Nevertheless, it is striking that IDavid et al.| 



J2009 



[20T0I ) 

find exactly this type of configuration in NGC 5044: few big 
cavities and many small 'weather-driven' bubbles, which 
appear to be radio quiet (i.e. not directly inflated by the jet). 

Overall, we conclude that the AGN outflows dynam- 
ics is particularly complex, as expected from the previous 
study in galaxy clusters. The main stages of evolution for 
the best models are: (i) a gentle continuous outflow, which 
sustains the cooling flow for the majority of time (or always 
in Be models) and which is barely observable in SBx maps; 
(ii) a phase of large cavity inflation, with cold rims and high 
density/SBx contrast (only for non-continuous models); (iii) 
small periods of quiescence in which the turbulence and vor- 
ticity promotes mixing and tiny weather-driven 'bubbles'. 

The presence of the above described features (big 
and small cavities, asymmetry of Fe abundance, low Mach 
shocks, etc. ), is consistent with several recent X-ray observa- 
tions of heated systems, and demonstrates the fundamental 
role of AGN driven outflows in the evolution of groups (and 
clusters). 



5 DISCUSSION 

In this paper we have proposed 3D numerical simulations of 
self-regulated AGN outflows in a galaxy group environment. 
We have followed the long term evolution of the interaction 
between two main characters: radiative cooling and feedback 
heating. The aim of this work is to understand the general 
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features of the feedback process to solve the coohng flow 
problem, i.e. quenching cooling rates, but, at the same time, 
preserving the global cool core structure. 

We tested two relevant models, which we had found to 
be successful in galaxy clusters ( [Gaspari et al.|[2011| |: cold 
feedback and Bondi feedback mechanisms. Galaxy groups, 
with masses ~ two order of magnitude less than big clusters 
are more fragile systems and the feedback must be conse- 
quently more delicate. In the following we discuss and sum- 
marize the main results obtained, along with merits and 
flaws of the models. 



5.1 Hot continuous Bondi (Be) 

In the first set of simulations the commonly adopted Bondi 
prescription is applied to evaluate the accretion rate onto 
the black hole. Although the Bondi scenario is quite naive, 
it can be treated as the template model for a typical self- 
regulated feedback directly linked to entropy. The efficiency, 
as free parameter, might just incorporate the discrepancies 
between the ideal Bondi model and the real dynamics. 

Overall, we found (similarly to clusters) that the Bondi 
accretion models are able to satisfy our strict requirements 
for a plausible evolution, with a mechanical efficiency in the 
range 5 x 10~^ — 10~^. As expected, the accretion rates are 
always two to three orders of magnitude lower than Edding- 
ton limit. This implies that the outflow 'touch' is always 
delicate and non-explosive, with power of the order of 10*'' 
erg s~' and velocity 10* km s~^. These values are consis- 
tent with numerous observations of AGN outflows through 



absorption lines (e.g. ,Nesvadba et al.||20lT K 

In the best models cooling rates are reduced below a few 
or fraction of M© yr~^, with a monotonic decrease in time. 
The density radial profiles of these models do not deviate too 
much from the observed ones. The same can be said of the 
temperature gradient, with a small overheating in the inner 
4 kpc nucleus. The last feature is probably associated with 
the continuous nature of the feedback: the Bondi formula is 
in fact based on the accretion of hot gas, whose reservoir 
never drains. 

A problem for this model may be the continuous injec- 
tion of energy, which carves a tunnel in the IGM. However, 
as shown in Section 4, this channel is usually very narrow 
and almost disappears in the surface brightness image (on 
the contrary of events associated with big bubbles). The 
important result is that a feeble continuous heating could 
operate, without being in contradiction with observations. 

The creation of cavities of tens kpc is another issue. 
However, we have noted that, during the increasing gen- 
eration of turbulence (due to the AGN), the jet becomes 
gradually more disturbed: eventually it will fragment, pro- 
ducing a series of small buoyant bubbles. This phenomenon 
is certainly enhanced by a cosmological evolution through 



merging and other large scale motions (e.g., Morsony et al 



20101 



5.2 Cold feedback (C) 

In galaxy clusters the triggering mechanism linked to the in- 
stantaneous AAfcoid induced a very powerful outburst, often 
super-Eddington. In much less massive galaxy groups such 
a powerful feedback would easily destroy the thermal struc- 
ture observed, producing negative temperature gradients up 
to 100 kpc. Therefore, successful models adopt efficiencies 
about an order of magnitude lower with respect to those 
used for clusters. With between 5 x 10~* and 10~^ cooling 
rates are still well below 10 per cent of the pure CP model, 
but now with consistent central densities and temperatures. 

The model is naturally impulsive, but this time the duty 
cycle is very high compared to the cluster evolution shown 
in Gil. In fact, for the majority of time the outflow is con- 
tinuous with moderate power, similar to the previous Bondi 
models. The AGN stays in a quiescent phase only 15-20 per 
cent of the time, after which it will ignites with outbursts up 
to 10*^ erg s~^. These powerful jets are essential for creating 
big bubbles of 10-20 kpc diameter. 

One of such event has been analysed in Section 4 (Fig. 
[6|. The bubble does not show high temperature within (e.g. 
10® K or more), like in clusters (Gil): it is very under- 
dense with respect to the ambient medium, with a relatively 
large contrast in the SBx map (> 50%). It is striking that 
the metal abundance is often highly asymmetric, enhanced 
along the outflow direction. All the previous features are 
commonly present in deep X-ray observations (McNamaral 



& Nulscn 2007] 


Baldi et al.| 


2009| GastaldoUo et al. 2009 


Kirkpatrick et a 




2009 


2011 


Rasmussen & Ponman(|2009 


David et al.||2009 


2010 





We have also checked that slightly different initial con- 
ditions, like an initial isothermal temperature profile and a 
NEW dark matter profile, do not alter the behaviour of the 
feedback, producing similar results. 



Overall, it is surprising that, in galaxy groups, models 
with cold feedback resemble, most of the time, simulations 
with Bondi accretion, with the merit of adding important 
features like X-ray cavities and ripples at larger radii (essen- 
tially weak shocks with Mach around 1.1 — 1.3). This fact 
clearly indicates the necessary properties of the feedback 
process, if outflows are the main mechanism: very frequent 
or almost continuous and gentle. 



5.3 Bondi with cold timing (Bi) 

Simulations using the third class of models, Bi, have been 
carried out to test the combined effect of Bondi accretion 
and cold triggering mode, in order to avoid Bondi continu- 
ous injection (intrinsic in the ideal case). In some prelimi- 
nary tests (here not shown), we found that a simple super- 
position of the two feedback schemes leads to unacceptable 
models, because the outcome is just a more powerful jet, 
which carves a larger and deeper continuous channel. 

Therefore, we tried another way to combine both meth- 
ods. Keeping Bondi accretion prototype (entropy-regulated) 
for the value of power, we assumed that when the gas begins 
to cool at the center (i.e. the inflow initiates), the feedback is 
activated. This will prevent too peaked instantaneous out- 
flow powers and, at the same time, the AGN activity ac- 
quires a sort of duty cycle. 

The results are very encouraging. The required Bondi 
efficiency for best models stays again in the range 5 x 10~^ — 
10~^. Assuming es < lO"'^, the evolution manifests the typ- 
ical flaws of a pure CF run. On the contrary, if eb > 10~^ the 
feedback becomes very explosive and easily 'burns' the cool 
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core. The best Bi models show instead radial profiles with 
gradients not far detached from the initial observed state, 
along with very low cooling rates (usually one or less solar 
masses per year). 

The look of the Bi simulations has been analysed, in- 
depth, in Section 4. The general features are almost identical 
to the previous discussed models. Extensive periods of a 
continuous Bondi-like injection (Pj ~ 10** erg s~^, wj ~ 
10* km s~^) produce a narrow heating channel, alternated 
with dormant AGN phases (~ 15%). As in the cold feedback 
computations, the short inactivity promotes the subsequent 
generation of outbursts with greater power (3-4 times more 
than the continuous phase), producing important observed 
features, such as tens of kpc cavities and oscillations in the 
flow quantities. 

It is relevant to note that, as in previous good cold mod- 
els, only a fraction of the injected mechanical energy is used 
to inflate cavities (usually under 50 per cenlQ). One of the 
main results, also pointed out in Gil for clusters, is that the 
success of a feedback crucially depends on how the energy is 
transferred to the surrounding medium, not only on the to- 
tal amount of available energy. Furthermore, the heated flow 
has a very complex dynamics, with numerous facets. After 
the mechanical energy is injected in the system, a part of 
energy inflates the bubble, but the rest generates the co- 
coon and weak shocks. When the jet is more continuous, 
the moderate energy released at intermediate radii induces a 
turbulence growing in time. These chaotic motions permits 
a more efficient circulation of the hot gas, promoting the 
deposition of energy through mixing at the base of the jet, 
where the cooling tends to dominate. Furthermore, while the 
system becomes more turbulent, the feeble jet can be frag- 
mented and can cyclically produce tiny 'bubbles', few kpc in 
size (albeit very faint in the SBx map). Their buoyant mo- 
tion may promote additional mixing and heating, helping 
the global feedback process. Quite simple feedback models 
still generate an amazing variety of heating processes. 

As a concluding remark for these models, the efficiency 
comparison underlines again the strict relationship with the 
cold best models, exactly around Ec ~ 5 x 10~* — 10""^. We 
can therefore affirm that the feedback in galaxy groups, able 
to quench cooling flow but at the same time preserve the cool 
core structure, is globally only of one type: moderate, gentle, 
often continuous and turbulent, with very short periods of 
inactivity and following outbursts. 



5.4 Rejected models 

It is also interesting to highlight the main features of unsuc- 
cessful models, in order to understand why some kinds of 
feedback do not work. 

First of all, while a small part of the total energy in- 
jected could be in thermal form, above a threshold fraction 
of ~ 50 % the center of the group becomes too hot with 
respect to the observations (see also Brighenti & Mathews 
2002, 2003). Injection of a modest amount of thermal energy 



^ E.g., the inflation energy of the bubble in Fig. [S] panel b, is 
APV ~ 4(2 X 10"" erg cm-^) (1.1 x 10^^ cm^) 9 X 10^'' erg, 
while injected mechanical energy is 2.45 X 10^* erg. Thus, inflation 
energy is just 37% of the total injected energy in that single event. 



can have the beneficial effect of disrupting the otherwise or- 
dered structure of the outflows (Section 3.5). This model 
intermittently eliminates the long tunnel, usually not ob- 
served in real systems. Moreover, the deposition of heat in 
the nucleus helps preventing the formation of a cold torus. 
Injection of thermal energy can be justified assuming that 
shocks heat the hot IGM in the very central region. 

Another rejected model considers the active outflowing 
mass directly linked to the real accreted mass onto BH (In- 
Out feedback. Section 3.6). The main problem here is that, 
even with very low /in (e.g. /out > 0.9), the jet velocity often 
exceeds 10^ km s~^. This unappealing result can in princi- 
ple be fixed adding a mass loading factor {rjfout', see also 



Cattaneo fc Teyssier 2007 Ostriker et al. 20101. However 



this assumption will mimic entrainment, just as our stan- 
dard outflow generation method described in Section 2.2. It 
is not surprising therefore that we recover acceptable results, 
similar to those described in Sections 5.1. 

Intermittent models (Section 3.4), with fixed jet power, 
are only partially acceptable, as was also found for galaxy 
clusters (Gil). The most successful run has relatively weak 
outflows (Pjct ~ 10~*PEdd, activated every 10 Myr, with 
a duration of 5 Myr each (model Int510m). We have also 
tested jets with Eddington power (~ 10*^ erg s~^). As ex- 
pected they will destroy the cool core generating systems 
very different from the real ones. 

Probably the major flaw of the intermittent models is 
their artificial nature. Cooling rates can be suppressed, com- 
pared to a pure CF model, but the temperature profiles 
oscillate strongly in the center. This further supports the 
common notion that heating must be self-regulated. 

Finally, we have briefly investigated the effect of jet pre- 
cession (not shown). Important parameters like inclination 
angles and revolution time of the jets are currently poorly 
known. Nevertheless, we found that the behaviour of cold 
gas at the base of the outflow, i.e. the formation of a small 
cold torus, is not greatly different from non-precessing mod- 
els. In fact, the heating channel will be again generated after 
the first episodes, getting a tilted orientation in respect to a 
normal z-axis cylindrical simulation. 

5.5 Comparison with galaxy clusters 

As expected, our simulations show that AGN heating has 
a deeper impact in galaxy groups than in clusters, because 
the same mechanical power per particle has a greater effect 
for a lighter, less bounded systems. 

As we have seen in the previous Sections, a jet with 
Eddington power can easily erase the entire thermody- 
namic structure of the cool core. The same applies for 
all the self-regulated simulations with very high efficiencies 
(> 5 X 10"^ for cold models; > 10"^ for Bondi). In Gil 
we showed that consistent cold models commonly produced 
super-Eddington outbursts (> 10*^ erg s~^), with a very low 
frequency (~ 10 per cent). Despite the huge power injected, 
the shock-heating phase could not destroy the observed den- 
sity and temperature gradients, and the cluster quickly re- 
stored the usual cool core. The only marked consequences 
were more ripples at larger radii and cavities, which in the 
early phase presented high internal energy and shocked rims. 

On the contrary, galaxy groups, with much less binding 
energy (per particle), can not recover from the same strong 
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impulsive feedback. In order to prevent such 'heating catas- 
trophe', cold feedback models in galaxy groups require an 
efficiency at least 5—10 times lower. The outcome is strik- 
ing, because the evolution resembles that of Bondi models, 
with low power jets (10** — 10*^ erg s~^) activated in a 
quasi-continuous gentle way. These outbursts occurring af- 
ter short quiescence periods (tens Myr) generate ~ 10 kpc 
cavities surrounded by relatively cold rims in the region close 
to the center of the group, similar to some observations (e.g., 
Gastaldello et al. 20091. Unlike in galaxy clusters, cavities 



tions of AGN outflows ( Crenshaw et al.]|2003 Morganti et 
al.|2005|[2007l|Nesvadba et al.|2008||2011| ). 



in groups present the usual low density contrast, with tem- 
peratures similar to the surrounding medium. 

The best models assuming continuous Bondi accretion, 
on the other hand, are equivalent in both clusters and groups 
(except for the efficiency, lower in groups). Accretion rates 
always follow a sub-Eddington regime. The moderate out- 
flows carve a narrow tunnel in the intergalactic medium; we 
showed, in this paper, that it is very faint in the SBx map 
and thus hard to resolve. The absence of inflated bubbles is 
still a riddle of this type of feedback. Notice that in galaxy 
clusters the low-power outflows (eb ~ 0.1) had more dif- 
ficulty in halting the massive cooling flow (tens M© yr~^), 
showing more peaked density profiles. Therefore, in clusters, 
cold accretion was slightly favoured as an acceptable model. 

Another common feature in both clusters and groups is 
the asymmetrical distribution of metals (mostly iron) pro- 
duced by the SNIa in the central galaxy. They tend to ac- 
cumulate along the jet axis, in qualitative agreement with 
recent observations (see Sec. 4). 

In all the consistent group models the sub-Eddington 
outflow powers imply that the total injected energy (a few 
10®^ erg for the full-space system) is always an order below 
the total 'available' BH energy. In fact, during the simulated 
evolution, the black hole mass increases at most 10 per cent 
in Bondi models. In the cold feedback model the lower ef- 
ficiencies imply a much higher AMbh, similarly to galaxy 
clusters. As noted in Gil, however, the real efficiency (still 
unknown) could just be higher, with a lower mass actually 
falling onto the BH. 

Our simple approximation of the real sub-pc accretion 
does not allow to answer why the efficiency in groups should 
be a factor 5-10 smaller than in clusters, aside the observa- 
tional requirement. It would be interesting to investigate this 
topic with a dedicated analysis. 



iii) The main feature of the continuous phase is a 
narrow tunnel (large a few kpc, up to 30 kpc long), which 
is almost undetectable in the SBx map. The tunnel is 
periodically fragmented by the increasing AGN-driven 
turbulence, especially when the BH is in a phase of low 
accretion. The products are tiny, kpc-size buoyant bubbles, 
which help the mixing of the IGM. These faint features are 
very difficult to detect, but seen in a few systems, most no- 
tably NGC 5044 ( |Buote et al.]2003l|Da^ et al.|2009|[2010 |. 



iv) Cold and hybrid accretion models show short 
quiescent periods (summing to a ~ 15% of the evolutionary 
time). Each quiescent period is followed by a relatively 
strong outburst, with power about an order of magnitude 
larger than the typical value during the continuous phase. 
The consequence is the creation of relatively large cavities 
(10 — 20 kpc), with high density contrast, in approximate 
pressure equilibrium with the ambient. The emission- 
weighted temperature in the cavity region, however, is 
not far above those of the surrounding medium. Feeble 
rims, relatively cold in the region near the center and 
slightly hotter in the more distant part, with somewhat 
higher iron abundance, are also present. These features are 



again reminescent of real systems (McNamara & Nulsen 
[2507| for a review, [Gastaldello et al.|2009[ [David et al.|2009 

v) The asymmetrical transport of metals, along jet-axis, 
up to 40 kpc, is clearly visible during the active phases. In 
quiescent periods, the moderate turbulence promotes metal 
diffusion. Recent observational data confirms this behaviour 
( Kirkpatrick et al.|20 09, 2011"] [David et al.|20ld l. 



vi) Finally, the global point of all our computations: 
the only possible way to heat a galaxy group, in order to 
suppress the grip of the cooling flow, appears to be through 
a quasi-continuous 'delicate touch', while in massive clus- 
ters rare and explosive strokes are mildly favoured over the 
weaker Bondi regulation. This quite different feedback evo- 
lution might explain the common observed discrepancies be- 
tween the two cosmological actors. Galaxy groups are indeed 
not purely scaled-down versions of galaxy clusters. 



6 CONCLUSIONS 

Here we summarize the main conclusions of our study. 

i) Feedback triggered by both Bondi accretion and 
cold mode accretion (or a combination) lead to successful 
self-regulated models. These schemes are able to quench 
cooling flow for many Gyr, preserving the observed 
thermal structure of relaxed groups. The required mechan- 
ical efficiencies for cold accretion models are in the range 
5 X 10"* - 10"^, while 5 x 10"^ - 10"^ for the Bondi scenario. 

ii) In the galaxy group, the two main feedback schemes 
generate similar flows. The global evolution is dominated 
by almost continuous outflows, with sub-Eddington power 
(10** — 10*^ erg s"'^) and velocities around 10* km s~^. 
These values are consistent with those based on observa- 
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